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ABSTRACT 

We present a sample of 2865 emission line galaxies with strong nebular He n A4686 
emissions in Sloan Digital Sky Survey Data Release 7 and use this sample to in- 
vestigate the origin of this line in star-forming galaxies. We show that star-forming 
galaxies and galaxies dominated by an active galactic nucleus form clearly separated 
branches in the Hen A4686/H/3 versus [Nil] A6584/Ha diagnostic diagram and derive 
an empirical classification scheme which separates the two classes. We also present an 
analysis of the physical properties of 189 star forming galaxies with strong HeiiA4686 
emissions. These star-forming galaxies provide constraints on the hard ionizing con- 
tinuum of massive stars. To make a quantitative comparison with observation we use 
photoionization models and examine how different stellar population models affect 
the predicted He II A4686 emission. We confirm previous findings that the models can 
predict He II A4686 emission only for instantaneous bursts of 20% solar metallicity or 
higher, and only for ages of ~ 4 — 5 Myr, the period when the extreme-ultraviolet 
continuum is dominated by emission from Wolf-Rayet stars. We find however that 83 
of the star-forming galaxies (40%) in our sample do not have Wolf-Rayet features in 
their spectra despite showing strong nebular He II A4686 emission. We discuss possible 
reasons for this and possible mechanisms for the He n A4686 emission in these galaxies. 
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1 INTRODUCTION 

The ionizing continuum of stars at A < 912 A is of ma- 
jor importance for interpreting emission line observations of 
galaxies because many of the strong lines observed in the 
spectra of galaxies, such as [O ill] A5007, [Ne ill] A3869 and 
Hell A4686, have ionization potentials in excess of 13.6 eV. 
Despite this importance we are severely limited by interstel- 
lar absorption in obs erving stellar spec tra in this spectral 
window directly (e.g. iHoare et alj fl993h Although we can 
get direct information at slightly lo nger wavelengths wit h 
space-based UV spectroscopy (e.g. ICrowther et al.l [2002). 
most of our knowledge about the A < 912A region is based 
on indirect evidence, even for solar metallicity. 

A promising way to indirectly obtain information on the 
stellar ionizing continuum is to compare emission line prop- 
erties (e.g. flux, equivalent width) to predictions from p ho- 
toionization code s such as CLOUDY l|Ferland et al.lfl998h or 
MAPPINGS III (| Allen et alj|2008l ). In practice these kinds 
of studies pr ovide modest constrai nts on stellar atmosphere 
models fe.g. ICrowther et al-l ^gg). However where predic- 



tions of models differ significantly, this approach can yield 
useful information. This is the approach we will adopt in 
this work, where we will make use of the HellA4686 nebular 
emission line to place constraints on stellar models and in 
particular on the ionization mechanism for this line. 

The presence of a nebular Hell A4686 line in the in- 
tegrated spectrum of a galaxy indicates the existence of 
sources of hard ionizing radiation as the ionization energy 
for He + is 54.4 eV (A w 228 A). This hard radiation can of 
course be produced by an active galactic nucleus (AGN), and 
most sources with luminous He II emission, in a flux limited 
sample, are indeed galaxies with an However the re- 

quired hard radiation can also be provided by stellar sources 
and He II A4686 emission is frequently seen in H n-galaxies. 
The line appears t o be associated with yo ung stellar popula- 
tions; for instance, iBergeron et al.l (Il997l ) proposed Of stars 
as the sources of HellA4686 emission in dwarf galaxies. Sub- 
sequent discussion has mostly focused on Wolf-Rayet (WR) 
stars, although the distinction between th e se two classes is 
rather blurred (e.g. iGrafener et al.ll201ll h ISchaererl (|l99(l 
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1 We will here not distinguish between the host galaxy and its 
nuclear power source so will refer to these galaxies as AGNs 
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see also lSchaerer fc Vaccal |l99Sl) ) showed that the hard ra- 
diation field of WR stars could provide a Roo d explanation 
of the nebular He n A4686 seen in H 1 1- galaxies. I Guseva et al.l 
(2000) tried to test this in a careful study of Hn-galaxies 
with prominent WR features. They were however unable to 
find WR features in 12 out of the 30 galaxies with nebular 
HenA4686 emission. The same lack of WR features in metal 
po or Blue Compact Dwar f (BCD) galaxies was pointed out 
bv lThuan fc Izotovl (|2005l ). They proposed that fast radia- 
tive shocks could b e responsible for this emission (see also 
iGarnett et al. 199ll ). 

Similar results were reported by iBrinchmann et all 

(2008, hereafter B08), who analysed a sample of strong 
emission line galaxies in the Sloan Digital Sky Survey 
jSDSS. York et all |2000| ) with He n A4686 emission. They 
showed that at least at metallicities of 12 + logO/H > 8, 
there appeared to be a close correlation between WR fea- 
tures in galaxies and the presence of Hell A4686 emission, 
but this appeared not to be so clear-cut at lower metallici- 
ties. 

This apparent lack of connection of Hell emission 
with the hard UV radiation from WR stars has also been 
seen in spatially resolved spectroscopic studies. iKehrig et. all 
(2008) performed an integral field spectroscopy study for 
the Hn galaxy IIZw70 and found that the region asso- 
ciated with nebular He 11 A4686 emission was a few arcsec 
offset from the region with detec ted WR features. More re - 
cently IKehrig et. all i|201ll ) and iNeugent fc Massevl (|201ll ) 
have presented studies of Hell A4686 emission in M33. 
Both studies find some regions with nebular Hell A4686 
emission that are not associated with WR stars (see also 
lHadfield fc Crowther 20071. lLopez-Sanchez fc Esteban 2010l 



and iMonreal-Ibero et al. 201Gi ). 

Thus a series of studies have shown that while 
He 11 A4686 emission frequently is found in association with 
WR stars, it appears not to be so in all cases, particu- 
larly at low metallicity. As mentioned above, possible addi- 
tional sources of high energy photons c ould be X-ray binaries 
(IGarnett et al. 199ll ). strong shocks (iDopita fc Sutherland! 
1996), low- level AGN activity and al ternative models for 
stellar evolution (|Yoon fc Langer|[2005l ). However the exist- 
ing studies do not show clear trends that allow us to distin- 
guish between these scenarios. 

Crucially the samples in most of the previous studies 
have not been selected specifically to study Hell emission 
lines. To make progress in understanding this puzzle it is 
important to have as large as possible sample of He 11 emit- 
ting galaxies to allow one to study the relationship between 
Hell emission and other physical properties. To this end 
we present here an analysis of emission line galaxies with 
strong He 11 A4686 em ission in SDSS Data Release 7 (DR7, 
lAbazaiian et al.ll2009l ). 

In section [2] we discuss the sample selection and care- 
fully account for AGN emission. The physical properties of 
the He 11 emitting galaxies are discussed in section [3] and 
the observed Hell A4686/H/? ratios are compared to model 
predictions in section[3] In section[S]we test these model pre- 
dictions and investigate whether the presence of He 11 A4686 
is associated with WR features. We find that low metallicity 
systems frequently do not show signs of WR stars. We dis- 
cuss possible explanations for this finding in section [6] and 
conclude in section [7] 



Measured line 


Number of spectra 


Fraction 


[On] A3726.3729 


243977 


16.51% 


[Nem] A3869 


39143 


2.65% 


H8 


92230 


6.24% 


[Nem] A3967 


18367 


1.24% 


He 


64120 


4.34% 


He 1 A4026 


6368 


0.43% 


[S 11] A4069 


1503 


0.10% 


H<5 


163534 


11.07% 


H7 


295528 


20.00% 


[O in] A4363 


7886 


0.53% 


He 1 A4472 


15563 


1.05% 


[Fe m] A4658 


2626 


0.18% 


He 11 A4685 


4034 


0.27% 


[Anv] A4711 


2893 


0.20% 


[Ar iv] A4740 


1074 


0.07% 


H/3 


458324 


31.02% 


[O in] A4959 


147734 


10.00% 


[O in] A5007 


339212 


22.96% 


[Ni] A5197 


1318 


0.09% 


[N 1] A5200 


758 


0.05% 


[CI in] A5518 


160 


0.01% 


[CI in] A5538 


190 


0.01% 


[Nil] A5755 


4888 


0.33% 


He 1 A5876 


75256 


5.09% 


[0 1] A6300 


119328 


8.08% 


[S in] A6312 


5552 


0.38% 


[U 1] A6363 


16100 


1.09/0 


[N 11] A6548 


361641 


24.48% 


Ha 


613338 


41.51% 


[Nil] A6584 


562659 


38.08% 


He 1 A6678 


19084 


1.29% 


[S 11] A6717 


410903 


27.81% 


[Sir] A6731 


341386 


23.11% 


He 1 A7065 


4232 


0.29% 


[Arm] A7135 


24034 


1.63% 


[0 11] A7318, 19,29,30 


659 


0.04% 



Table 1 . The table shows number of spectra that have the indi- 
cated line detected at S/N > 5.5. The total number of analysed 
spectra is 1,477,411. 



2 DATA 

Our sample is based on galaxy spectra from SDSS DR7 
which cover a wavelength range of 3800-9200 A. The spec- 
tra were analysed using the methodology discussed in 
iTremonti et al. 20041 (see also Brinchmann et al. 2004) to 
provide accurate continuum subtraction. All emission line 
sources were additionally analysed using the pipeline dis- 
cussed in B08 to measure a wider gamut of emission lines. 
For each galaxy we measure 40 lines, These lines and the 
number of spectra that show these lines with S/N > 5.5 are 
summarised in Table [T] 

Our concern in this paper is not to analyse a volume- 
or magnitude- limited sample of galaxies, we therefore do not 
impose a redshift cut nor a magnitude limit. Since the blue 
wavelength cut-off of the SDSS spectrograph is ~ 3800A, the 
[0 11] A3727, 3729 doublet falls outside the spectral range for 
z < 0.02. This is a concern, because as we will see later 55% 
of our final sample fall in this region of redshift space and 
their oxygen abundances are therefore somewhat uncertain. 
When possible we use the [O n] A7318 — 7330 quadruplet 
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instead (|Kniazev et al.| [2004). but as this is a fairly weak 
line and falls in a region with significant sky emission, we 
cannot always make use of this line. 

2.1 Sample selection and classification 

We select our sample requiring signal to noise ratio > 5.5 
in Hen A4686, the resulting data set is given in Table [2]Q. 
When the width of the He II line is consistent with that of 
the strong forbidden lines, we make the assumption that it 
has a nebular origin. Given that Hell lines from individual 
WR stars typically are considerably broader than the for- 
bidden lines in galaxies (e.g. B08), we feel this is a reason- 
able assumption. In addition we require a S/N> 3 in each 
of H/3, [O ill] A5007, Ha and [N n] A6584 emission lines to 
reliably classify our galaxies (Brinchmann et al. 2004, here- 
after B04) . The resulting sample contains 2865 spectra with 
strong nebular HellA4686 emission. In parallel, the spectra 
of the sample galaxies are examined for WR signatures us- 
ing the approach discussed by B08. This resulted in a total 
sample of 385 spectra with likely and secure WR features 
(Class 1, 2 and 3 from B08). While we do not discuss the 
sample of all WR galaxies in DR7 with Class 1-3 in detail 
here, we note that it intersects that of the He II sample but 
is not a strict subset of it (see Table [3}. 

Figure [JJ shows the redshift distribution for the frac- 
tion of the Hell sample in the SDSS as a shaded grey his- 
togram. The cut-off at z ~ 0.4 is due to Ha falling outside 
the spectrograph range. The red histogram shows the red- 
shift distribution for the fraction of just the star-forming 
(SF) galaxies in the SDSS showing HellA4686 emission (see 
below for a discussion of the classification). For each class we 
have divided the number of He II emitting galaxies in that 
redshift bin by the number of similarly classified galaxies in 
the parent sample (SDSS galaxies that have S/N> 3 in each 
of Hj8, [Om]A5007, Ha and [Nll]A6584) in that redshift bin. 
A constant value would therefore indicate a similar redshift 
distribution of the Hell sample and the parent sample. It 
is clear from this that full Hell sample closely follows the 
overall distribution of the SDSS, but that the star-forming 
galaxies with He II A4686 emission are predominantly found 
at low redshift. We can also see that less than 2% of all 
galaxies, and less than 1% of the star-forming galaxies in 
the SDSS DR7 show Hen A4686 emission in their spectra. 

To classify the dominant ionization source in 
each galaxy we fo llow previous studies in using the 
iBaldwin et al.l |l98ll . BPT) line ratio diagnostic diagram 
of [Om] A5007/H/3 versus [Nil] A6584/Ha as our starting 
point (Figure ED). As has been discussed extensively (e.g 



0.020 



Terlevich et~aITll99ll; iKewlev et all 12001. hereafter KeOl 



Kauffmann et al. 2003, hereafter Ka03; IKewlev et al. 2006 



Stasinska et al.l 20061 ') this diagram allows a separation of 
AGN and star-forming galaxies because of their signifi- 
cantly different ionizing spectra, typically leading to high 
[Om] A5007/H/3 and [Nil] A6584/Ha when an AGN is 
dominating the output of ionizing photons. KeOl used a 
combination of stellar population synthesis and photoion- 
ization models to compute a theoretical maximum starburst 



2 The full table of 3292 spectra is available in electronic form in 
http: / /www. strw.leidenuniv.nl/~shirazi/SB011/ 




0.2 0.3 
Redshift 



Figure 1. The redshift distribution of the full Hell sample, rel- 
ative to that of the parent sample is shown as a shaded grey 
histogram. The red histogram shows the redshift distribution of 
the He II SF sample relative to that of all star- forming galaxies in 
the parent sample; the parent sample consisting of those galaxies 
in the SDSS DR7 having H/3, [Om] A5007, Ha and [Nil] A6584 
detected at S/N> 3. In comparison with the SDSS, it is clear 
that the galaxies in our SF sample have a redshift distribution 
strongly shifted toward low redshift. 



line that isolates objects whose emission line ratios can be 
accounted for by photoionization by massive stars (below 
and to the left of the curve) from those where some other 
source of ionization is required. Ka03 defined an empirical 
upper limit to the Hn region sequence of SDSS galaxies 
in the BPT diagram. The region lying between these two 
lines represents objects more naturally explained as having 
a composite spectrum combining H II region emission 
with a harder ionizing source. As we will see later, this 
interpretation is further corroborated by the He IIA4686/H/3 
ratios we find for our sample. 

In the present study, we adopt a two-stage classifica- 
tion methodology. We start out by classifying all galaxies 
using the BPT diagram, and we will then refine our classifi- 
cations for a subset of the galaxies using detailed inspection 
of the spectra and the Hell A4686 line properties. For the 
initial classification we adopt a similar methodology to that 
of B04 and use the separation criteria defined by KeOl and 
Ka03 to divide galaxies into different classes. Galaxies which 
are distributed above the KeOl dividing line are considered 
AGNs, Galaxies between the KeOl and Ka03 limits have 
composite classification, which means their source of ioniz- 
ing radiation could be a combination of star formation and 
AGN activity. Finally, galaxies below the Ka03 line have a 
SF classification, which as we will see might be modified 
subsequently. 

Figure [5] shows the BPT diagram for our sample, the 
KeOl and Ka03 classification lines are shown as solid and 
dotted lines, respectively. The distribution of all emission 
line galaxies in the SDSS with S/N > 3 in H/3, [O ill] A5007, 
Ha and [N n] A6584 is shown as a grey-scale 2D distribu- 
tion, where the grey-scale shows the logarithm of the num- 
ber of galaxies in each bin. Blue circles show star-forming 
galaxies, triangles show composite galaxies and red squares 
mark AGNs. The grey dashed-dotted line shows the N2 
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1.5 




Figure 2. This plot shows the BPT diagnostic diagram for the 
sample. The Ka03 classification line is shown as a dotted line and 
KeOl classification line as a solid line. The grey dashcd-dotted line 
shows the N2 PP04 metallicity calibration for 12 + log O/H = 8.2. 
The distribution of emission line galaxies in the SDSS is shown 
by the gray-scale 2D distribution where the grey-scale shows the 
logarithm of the number of galaxies in each bin. Blue circles show 
star-forming galaxies, triangles show composite galaxies and red 
squares mark AGNs. As discussed in the text, for some galaxies 
the classification has been adjusted which is why some objects 
in the star-forming region are classified as AGN. We mark these 
galaxies with a blue plus over the red square. 



l|Pettini fc Page! 2004 hereafter PP04) metallicity calibra- 
tion for 12 + log O/H = 8.2 ([N n] A6584/Ha = -1.2). These 
classifications are the final ones and incorporate further in- 
formation as described in the following. 

While the BPT diagram is a useful classification di- 
agram, it is not particularly sensitive to low levels of 
AGN contamination and some progress can be made by 
including lines orig inating in the mostly neutral ISM (e.g. 
iKewlev et alll2006h . For our purposes we however need to 
be very confident in the lack of AGNs in our sample and we 
will use the Hell A4686/H/3 ratio for this purpose. 

As remarked earlier, only photons with energy in excess 
of 54.4 eV can ionize He + and thence produce the He IIA4686 
recombination line. If we consider single stella r population 
models from Starburst99 (|Leitherer et al.lll999h at an age of 
2 Myr (so that all stars are on the main sequence), we find 
that the integrated spectrum of this population typically 
contain 2-3 orders of magnitude fewer photons at this energy 
than at the energy required to ionize + (35.5 eV), which 
is needed to produce the [O ill] A5007 line from collisionally 
excited +2 . This line is therefore a very sensitive probe of 
AGN activity, particularly if used in conjunction with other 
line ratios. We therefore make use of the Hell A4686/H/3 vs 
[Nil] A6584/Ha diagram to further refine the classification 
of our sources. We show this diagram in Figure [3] where 
the symbols and colours are the same as in Figure [2] Note 
that there is generally a very significant offset between the 
star-forming galaxies and AGNs in this diagram, in contrast 
to the gradual transition in much of the BPT diagram. To 
make a quantitative separation, we draw a random set of 
star-forming galaxies from the SDSS and gradually add their 
emission line fluxes to that of a random set of AGNs. We 




Log [N II] X6584/Ha 

Figure 3. This plot shows our sample in the HenA4686/H/3 ver- 
sus [No] A6584/Ho diagnostic diagram. The dotted line shows 
an empirical line separating AGN and composite objects from 
star-forming galaxies. The grey dashed-dotted line shows the N2 
PP04 metallicity calibration for 12 + log O/H = 8.2. Symbols are 
the same as the BPT diagram. As He II A4686 have a higher ion- 
ization potential in comparison to [O III] A5007 and is much less 
sensitive to the electron temperature, we can clearly see separa- 
tion between the classes in this diagram. The lines with colour 
gradients in the figure are simulated fluxes drawn from adding a 
random set of star-forming galaxy from the SDSS with emission 
line fluxes of a random set of AGN spectrum. The colouring of 
the lines corresponds to the fraction of the spectrum contributed 
by the star-forming galaxy, 1 — /, as indicated by the colour bar 
on the side. The solid line shows the theoretical upper limit for 
Hen A4686/H/3 ratio. 

quantify this by the variable / which is defined to be the 
fraction of the total H/3 flux comes from the AGN. The total 
flux is therefore / x AGNfjaix + (1 — /) x SFfujx- Changing 
/ will trace out a path in the diagnostic diagram in Figure [3] 
as illustrated by the lines with colour gradients in the figure. 
The colouring of the lines corresponds to 1 — /, as indicated 
by the colour bar on the side. We repeat this process for a 
thousand AGN and SF objects located in different bins in 
the He II /H/3 diagram. 

Based on this analysis we find that there is a well defined 
locus where 10% percent of the He II A4686 flux comes from 
an AGN (in this case ~ 1% of the H/3 flux comes from the 
AGN, see section I2T21 for further details). A good fit to this 
relation is given by: 

log( FHe ° A4686 )=-1.22+ =r-± , (1) 

F W 8.92 log( [N ^ 6584 ) + 1.32 

which is shown as a dotted line in Figure [3] The solid line 
shows a theoretical maximum starburst line, similar to the 
KeOl line in the BPT diagram — we discuss this further 
below. 

Finally, we look at the spectra of galaxies that we would 
classify as star-forming on the basis of their location in 
the BPT diagram, but that are offset from the rest of the 
star-forming galaxies in the Hell /H/3 diagram and check 
whether they show AGN features such as broad Balmer 
lines, strong NevA3426, Fell emission or if they show a sim- 
ilar [O in] A4363/H7 ratio to that of AGNs. If some of these 
features are present, we change the classification from star- 
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Sample 


Hen 


Hen + WR 


Total 


2865 


385 


AGN 


2474 


234 


Star forming 


199 


116 


Composite 


179 


35 



Table 3. An overall numerical summary of the Hell sample. See 
scction[2]for a summary of the selection and classification details. 
See section Q>] for details on the WR classification. 



forming to AGN (these objects are plotted as red squares 
containing a blue cross in Figures [2] and [3} . This happens 
for 127 (39%) galaxies. This is a conservative approach as 
we would exclude star-forming galaxies with strong outflows 
and hence a broad base to the Balmer lines for instance. 

To estimate the maximum starburst line in Figure [3j 
we adopted the Chariot & Longhetti (2001, hereafter CL01) 
models. These combine evolving stellar populations models 
from Bruzual & Chariot (unpublished BC00 models using 
Padova 1994 tracks, see Bruzual & Chariot 2003 (BC03) for 
the current models) with the photoionization code Cloudy 
l|Ferland et alJll99S|) and adop t the s imple dust attenuation 
prescription of ICharlot fc Falll (|2000h . The main model pa- 
rameters for our calculations are the metallicity Z, the ion- 
ization parameter, U, the dust attenuation tv and the dust- 
to-metal ratio, £, the model parameters used are given in 
Table [4] see B04 for a more detailed discussion. We use the 
CL01 Single Stellar Population (SSP) models since these 
achieve the highest possible Hell A4686/H/3 values. We then 
identify the maximum ratio reached by the different models 
and use this upper envelope to define the maximum star- 
burst line shown as a solid line in Figure [3] 

This combination of classification methods means that 
there is not a one-to-one mapping between the location of 
an object in a diagnostic diagram and its final classification, 
as is clear from Figure [5] and We also mention that we do 
not change classes for galaxies classified as AGNs or com- 
posites using the BPT diagram, but which fall within the 
star-forming region in the He II /H/3 diagram and this is the 
reason why there are 7 AGN below our star-formation-AGN 
dividing line. 

By contrasting Figure [5] and [3] we can make a couple of 
interesting observations. The first is that while in the BPT 
diagram we see a steady increase in [O ill] A5007/H/3 with de- 
creasing [N n] A6584/Hq, in Figure[3]we see no major change 
in HenA4686/H/3 with [Nil] A6584/Ha. This indicates that 
He II A4686/H/3 and consequently the ionizing spectrum of 
stars at A < 228 A vary only weakly with metallicity. The 
second difference between the two plots is that the He II /H/3 
diagram can separate star-forming galaxies from composite 
galaxies better than BPT diagram since the He II A4686/H/3 
ratio is more sensitive to the hardness of the ionizing source 
than [O in] A5007/H/3. We can use this to further support 
our supposition that the gas in the galaxies falling between 
the Ka03 and KeOl lines in the BPT diagrams is ionized by 
a combination of stars and an AGN — while they are adja- 
cent to the star-forming sequence in the BPT diagram they 
nearly all clearly separate from the star-forming sequence 
in the Hell /H/3 diagram, corresponding to an AGN con- 




o.o t 

0.0 



0.4 0.6 
Fraction of H/S flux from AGN 



1.0 
0.9 

o.a 

0.7 

0.6 

0.5 
0.4 
0.3 
0.2 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 

Log [N II] A6584/Ha 

Figure 4. The upper panel shows the sensitivity of different lines 
to the presence of an AGN as a function of the fraction of the H/3 
flux originating from an AGN. The grey shading indicates the 
region where a typical galaxy would be classified as star-forming 
in the BPT diagram (see text for details). We see that we might 
classify galaxies as star-forming even if 10% of the H/3 flux comes 
from AGN. The bottom panel illustrates the method used and 
shows the path traced by a galaxy in the BPT diagram as an 
increasing amount of AGN light is added to a star-forming galaxy 
spectrum. The colouring of the line corresponds to the fraction 
of the H/3 line flux contributed by the star-forming galaxy, 1 — /, 
as indicated by the colour bar on the side. The dashed line is 
the Ka03 classification line, the intersection of this line with the 
trajectory of the simulated flux is at / ss 0.1. 



tribution to Hell A4686 > 50%. Thus referring to these as 
composite objects appear to be justified. 



2.2 AGN contamination estimation 

In view of the clear separation of the AGN and star- 
formation branches in Figure [3] and our sensitivity to 
low-level AGN contamination, it is beneficial to study 
the impact of an AGN on the line ratios in more de- 
tail. This has been di scussed in previous studies (e.g. B04, 
Stasihska et all l|2006t» but here we extend those efforts to 
include Hen A4686. 

We focus our attention on the BPT diagram as it is most 
widely used for emission line classification. We follow the 
same approach as in the previous section of adding gradually 
more of an AGN emission line spectrum with He II A4686 in 
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Parameter 



Range 



Z, The metallicity 
U, The ionization parameter 
Ty, The total dust attenuation 
£, The dust- to- metal ratio 



-1 < logZ/Z e < 0.6, 24 steps 
-4.0 < logll < -2.0, 33 steps 
0.01 <t v < 4.0, 24 steps 
0.1<{< 0.5, 9 steps 



Table 4. The model grid used for the present work. We calculate 
this both for a constant star formation history at t = 10 s yrs as 
well as for an SSP. 



emission at S/N> 3, to a star- forming one, and find where it 
intersects the Ka03 line (see bottom panel of Figure 4). At 
this point the galaxy would cease to be classified as a star- 
forming galaxy. We repeat this for a total of f 0,000 random 
combinations of spectra. The median AGN contribution to 
H/3 at the point where a galaxy ceases to be classified as 
star- forming is « 10%. 

Figure 4 shows the region where a median galaxy would 
be classified as star-forming as the gray shaded region. On 
top of this we show the median trend for the fraction of flux 
in the indicated line as a function of the fraction of the H/3 
flux coming from an AGN (shown for reference as the dashed 
diagonal line). We note that the exact shape and location of 
the [N n] A6584 and [O ill] A5007 lines does depend somewhat 
on the AGN sample chosen but the qualitative trend remains 
the same. 

What this figure shows, is firstly the well-known result 
that some lines are more sensitive to the presence of an 
AGN than others. As mentioned before, the Balmer lines 
are expected to have less than 10% contribution from an 
AGN, while the [Nil] A6584 line can have more than 30% of 
its flux coming from an AGN, putting in qu estion its use as 
an ab undance indicator on its own (see also [Stasirisk a~et al.l 
(2006)). But for our purposes, it is more important to note 
that by adopting a classification based on the BPT diagram, 
we would classify a galaxy as star-forming even when ~ 65% 
of its Hell A4686 flux would come from an AGN. 

We can now combine this with our previous result in 
Figure where we found that only 10% of the He II A4686 
emission comes from an AGN in our refined star-forming 
sample. Applying this to Figure [4] we conclude that less 
than 1% of the H/3 flux comes from an AGN and the other 
lines will also only have very small contributions from an 
AGN implying we have a quite pure star-forming sample. 

Our final sample of He II star-forming galaxies consists 
of 189 star-forming galaxies (199 spectra which have been 
summarised in Table [5]). As mentioned above we have also 
checked these spectra for the presence of WR signatures. We 
will return to a detailed discussion of this in section [S] but 
will use the result of this classification in the following plots. 



3 PHYSICAL PROPERTIES OF THE SAMPLE 

While the majority of the galaxies in our sample have some 
physical parameters in the MPA-JHU value added cata- 
logue^, our galaxies are sufficiently extreme that we need to 
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Figure 5. This plot shows Hell A4686/H/3 as a function of oxy- 
gen abundance (see text for details on the calculations). Oxygen 
abundances derived using the direct method are shown by black 
circles, we can not calculate oxygen abundances for 13 objects 
with this method as no [O II] lines are available for them. Oxygen 
abundances derived from fits to the CL01 models are shown by 
red circles while the blue squares show the 03N2 oxygen abun- 
dance estimates. Note that the qualitative trends are similar, but 
the 03N2 estimator does not reach as high O/H values as the 
other two models. At low metallicity the three methods are in 
good agreement. 



rederive some properties and add some physical parameters 
to what is in the MPA-JHU catalogue. 

For the calculation of physical parameters we will adopt 
the Bayesian methodology outlined by Ka03 and B04. For 
each model we calculate the probability of that model given 
the data assuming Gaussian noise and obtain the Probability 
Distribution Function (PDF) of every parameter of interest 
by marginalisation over all other parameters (see Appendix 
A). We take the median value of each PDF to be the best 
estimate of a given parameter. 

3.1 Mass measurements 

The standard SDSS pipeline often segment nearby actively 
star-forming galaxies incorrectly, thus we need to redo the 
photometry of our galaxies. We do this using the Graph- 
ical Astronomy and Image Analysis Tools (GAIAQ). Most 
of our galaxies have strong emission lines within some of 
the broad-band filters. Prior to fitting we therefore correct 
the magnitudes for emission line contributions by assum- 
ing that the relative contribution of the lines found in the 
SDSS fiber spectra is applicable to the galaxy as a whole. As 
most galaxies in our sample appear to have a uniformly blue 
colour, presumably due to active star formation, we expect 
this to be a reasonable assumption. 

Stellar masses are calculated as outlined above, by fit- 
ting a large grid of stochastic models to the SDSS u, g, r, i, z 
band photometry. The grid contains pre-calculated spectra 
for a set of 100,000 different star formation histories using 



i syr 

cepts of <|Gallazzi et ai1l2~005l . 120081 ) 



http : //www.mpa-garching .mpg. de/SDSS/DR7 



http : //astro . dur . ac .uk/7,7Epdraper/gaia/gaia .html 
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Figure 6. T he contours show the mass-metallicity relation for 
SDSS galaxies llTremonti et al. 20041) . The present sample of star- 
forming galaxies with WR features is shown by blue circles, while 
the red triangles show the locations of those that do not show 
WR features. At low mass we see that our sample is offset from 
the rest bulk of the of SDSS but overall they sample much the 
same region. 



3.2 Emission line derived parameters 

We use the CL01 model to analyse the emission lines in our 
sample. We adopt a constant star formation history (SFH) 
and use the same grid used by B04 (see Appendix A and 
B04 for further details). In total the model grid used for the 
fits have 2 x 10 5 different models. 

The main quantity of interest for the present discus- 
sion is the oxygen abundance, quantified as 12 + logO/H. 
As there are significant differe nces between methods fo r es- 
timating oxygen abundance dKewlev fc Ellison] 120081 ). we 
have complemented the estimate from the CL01 method 
with two independent methods: Firstly, we estimate gas- 
phase oxygen abundances with the empirically calibrated 
estimators proposed by PP04. They used the line ratios of 
[Om] A5007/H/3/[Nn] A6584/H«, the 03N2 method, and 
[Nil] A6584/Ha, the N2 method, as abundance indicators. 
For all objects with detected [O ill] A4363 with S/N > 3, we 
also use the T e metho d, or direct method , using the fitting 
formulae provided by llzotov et alj (|2006h to estimate the 
oxygen abundances. For those objects without [O ill] A4363, 
we adopt the electron temperature estimates from the CL01 
models fits for the direct method calculation. Whenever we 
do not have [O n] A3727, 3729, we use the [O n] A7320, 7330 
lines to calculate abundances. However, for 13 objects we 
are unable to use the direct method for estimating oxygen 
abundance as none of the [O n] lines are available. 

Figure [5] compares different abundance indicators in the 
Hell A4686/H/3 flux ratio versus oxygen abundance plane. 
The oxygen abundances derived using the direct method are 
shown by black circles while those derived from the fit to the 
CL01 model are shown by red circles, the blue squares show 
03N2 oxygen abundances. The main conclusion we can draw 
from this comparison is that all methods agree well at low 
metallicity while at high metallicity the trends are similar 
but the 03N2 estimator reaches a lower maximum O/H. 
For concreteness we will adopt the CL01 estimates for the 



Figure 7. This figure shows the logarithm of the ratio of the 
model prediction for Hell A4686/H/3 to the observed ratio as a 
function of oxygen abundance. Red triangles show the ratio for 
objects without WR features. It is clear that there is good agree- 
ment in the range 8.4 < 12 + log O/H < 8.8, but at lower metal- 
licity, the discrepancy between model and observations can be up 
to and order of magnitude. At higher than solar metallicity where 
an AGN contribution to the He II A4686 flux is more likely we see 
that model also fail to predict the same ratio as the observed 
value. 



remainder of the paper but as our main focus will be on the 
low metallicity region, our results are robust to the estimator 
chosen. 

Figure HJ] shows the mass-metallicity relation for the 
sample compared to the mass-metallicity relati on for all 
star-forming SDSS galaxies (|Tremonti et al. 20041 ) shown as 
a contour. The sample galaxies with WR features in the 
spectra are shown as blue circles, while those that do not 
show WR features are plotted as red triangles, we will use 
the same symbols in the following. Overall there is a reason- 
able agreement with the main SDSS sample except for an 
offset towards slightly lower metallicity at a fixed mass at 
low masses. 



4 MODEL PREDICTIONS 

In the previous section we carried out an empirical anal- 
ysis of the properties of star-forming galaxies in the SDSS 
which show strong nebular He II A4686 emission in their inte- 
grated spectra. Now we will build on the preceding to explore 
whether current stellar models can be used to explain the 
He IIA4686 emission seen in the spectra of these galaxies. We 
start by predicting nebular He II A4686 emission for galaxies 
in our sample with the CL01 model. Then we change the 
stellar population model and explore the effect of changing 
the model on the predicted He II A4686 line flux. 



4.1 CL01 predictions for nebular Hell emission 

We follow the same procedure as in the calculation of PDFs 
for the galaxy parameters in the previous section and cal- 
culate the likelihood of the model for each object in our 
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Figure 8. This figure shows the spectral energy distribution (SED) of an instantaneous burst calculated with Starburst99. Each panel 
corresponds to one metallicity and shows the SED for a range of ages as indicated. Note that the appearance of WR stars 4 Myr after 
the burst results in a much harder UV continuum. After 5-6 Myr the WR stars disappear and the UV continuum rapidly fades. Also 
note that in these models, the lowest metallicity SED does not show a significant WR phase. 



sample by fitting the CL01 grid of models to the five im- 
portant [O n] A3727, 3729, H/3, [O III] A5007, Ha, [N n] A6584 
emission lines. 

We now want to see whether the models that reproduce 
the main strong lines in the optical spectrum also repro- 
duce the Hell A4686 emission line strength. We build the 
likelihood distribution of He II A4686 flux for each galaxy in 
the same way as before by weighting the He II A4686 flux 
in each model by the probability of that model. We take 
the median of the likelihood distribution as a prediction for 
nebular He II A4686 emission and the associated confidence 
interval to be the 16th-84th percentile range. We follow the 
same approach to estimate the He II A4686/H/3 ratio. 

In Figure [7] we show the logarithm of the ratio of 
this model prediction for Hell A4686/H/3 to the observed 
Hell A4686/H/3 ratio as a function of oxygen abundance. 
It is clear that there is acceptable agreement between the 
model predictions and the observations in the range 8.4 < 
12+log O/H < 8.8, but a model that can reproduce most the 
strong lines in the spectrum well, predicts up to one order 
of magnitude lower Hell A4686/H/3 ratio than the observed 
value for some objects at lower metallicities. We also see a 
deviation at high metallicity — in this regime an AGN con- 



tribution to the Hell line flux is more likely, both because 
the galaxies are more massive, and also because the star 
formation- AGN separation is more gradual in this regime. 
We will not discuss this mismatch further here. 

The population synthesis model used in the CL01 mod- 
els approximate the WR emission as black bodies at their 
effective temperature, note in passing that this is not the 
case in the current BC03 models. This will overestimate 
the hardness of the ionizing spectra compared to models 
that consider more soph isticated WR atmosp here models 
such as e.g. Starburs t99 (|Leitherer et al.lll999h and BPASS 
l|Eldridge et al.|[2008l ). Given the relatively simple treatment 
of the WR phases in the CL01 models, one might be con- 
cerned that the failure to match the data is due to an in- 
herent weakness of the models. In the following we therefore 
look at the effect of different stellar evolution and atmo- 
sphere models on the prediction of He II A4686 emission line 
strengths. 



4.2 Starburst99 predictions 

To better understand the origin of the Hell A4686 emis- 
sion and its dependence on the stellar models adopted, we 
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Figure 9. Each panel compare the calculated SED of two SFHs with starburst99 for a different metallicity, solid lines show SEDs of 
instantaneous burst and dotted lines show SEDs with constant star formation rate and different colours show different burst ages. Fluxes 
have been normalised to the flux at 912 A. We can see although models with constant star formation form WR stars continuously after 
4 Myr, the overall shape of the UV continuum is softer than in the instantaneous burst models because the continuous formation of 
luminous O stars softens the extreme UV spectrum for a fixed rate of hydrogen ionizing photons. 



use the latest versi on (6.0) of the spectral sy nthesis code 
Starburst99 (Stb99, lLeItherer et al1ll999l . 120101 ) to calculate 
spectral energy distributions (SEDs) predictions for a range 
of ages and metal abundances. We calculate models with an 
instantaneous burst and a constant star formation history 
with a Kroupa IMF. We have explored a range of stellar 
evolution models but the differences are small so we only 
show the results of one model. For this we adopt the Padova 
AGB evolutio nary tracks comb ined with Pauldrach/Hillier 
atmospheres JSmith et al. 20021), St b99 uses O star model 
spectra from iPauldrach et al.l (l200ll) a nd W R model spec- 
tra from the code of iHillier fc Millerl (|l998l ). The models 
include stellar and nebular continuum. We create mod- 
els with different metallicities (0.0004, 0.004, 0.008, 0.02, 
0.05), where the referenc e solar metallicity is Zq = 0.02 
ijAnders fe Grevesse 19891 ) . We do not consider dust and run 
the models up to 100 Myr with time-steps of 10 years. Fig- 
ure [8] shows the resulting SEDs for an instantaneous burst 
with a range of metallicities. Each panel corresponds to one 
metallicity as indicated and shows the time of the SED for 
six burst ages. The plots show clearly that the appearance of 
WR stars, 4 Myr after the burst, results in a harder UV con- 



tinuum shortwards of the He + ionizing edge at 228 A. After 
~ 5 Myr the WR stars disappear and the UV continuum be- 
comes softer. This is in good agreement with the discussion 
in Schaerer & Vacca (1998, see their Figure 9), which is nat- 
ural as those models lie at core of the WR modeling in Stb99. 
Figure[9]compares the SEDs of an instantaneous burst (solid 
line) with that of a continuous star formation model (dotted 
line). The SEDs have been normalised at 912A, so have the 
same amount of hydrogen ionizing photons. The continuous 
star formation models form WR stars continuously after 3 
Myr, but the overall shape of the UV continuum is softer 
than for instantaneous burst models because of the continu- 
ous formation of luminous O stars which dilute the SED for 
a given total mass. 

To calculate the emission lines, we use the UV con- 
tinuum generated by the Stb99 models as an input to 
the photoionization code Cloudy (version c08. lFerland et all 
Il99ct ). For each time step, ionization bounded models are 
calculated by varying the ionization parameter log U = 
— 2, —3, —4, where for consistency with the CL01 mod- 
els we calculate the ionization parameter at the edge of 
the Stromgren sphere, and a constant hydrogen density 
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Figure 10. The left panel shows the starburst99 instantaneous burst model prediction for the He II A4686/H/3 ratio for different ionization 
parameters at a fixed age of 4 Myr. The triangles along each line correspond to the metallicities (0.02, 0.2, 0.4, 1) Zq and our sample 
galaxies are shown as purple circles. The right hand panel shows the same for a range of different metallicities and ages with logU set to 
—2. The triangles along each model line correspond to the burst ages (3, 4, 5, 6) Myr. Note that the lowest metallicity model is unable to 
cover the observational data. The grey dashed-dotted lines in each panel show the N2 PP04 metallicity calibration for 12 + log O/H = 8.2 
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Figure 11. The plots show the prediction of instantaneous burst binary model (BPASS code) for Hell A4686/H/3 ratio for different 
metallicities (0.05, 0.2, 0.4,1) zq and burst ages. The grey dashed-dotted line shows the N2 PP04 metallicity calibration for 12 + log O/H = 
8.2. 
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Figure 12. The plots show the prediction of instantaneous burst single-star model using the BPASS code for the Hen A4686/H/3 ratio 
for different metallicities, (0.05, 0.2, 0.4, 1) Zq, and burst ages. The grey dashed-dotted line shows the N2 PP04 metallicity calibration 
for 12 + logO/H = 8.2. 



of logriH/cm~ 3 = 2.5. Th is range in ionization p arame ter 
spans the range found by IStasihska fc Leithererl l|l996l ) in 
their analysis of intensely star-forming galaxies (approxi- 
mately —3.5 < logU < —2.5). Figure [TO] shows the calcu- 
lated HellA4686/H/3 ratio versus [Nil] A6584/Ha for differ- 
ent metallicities and different ionization parameters in the 
left panel and different burst ages in the right panel. The 
triangles along each model line correspond to the metallici- 
ties (0.02, 0.2, 0.4, 1) Zq in the left panel and to ages of 3, 
4, 5, and 6 Myr in the right panel. In the left panel we fix 
the the age to 4 Myr, and in the right we set log U = —2. 
The galaxies in our sample are shown as filled purple cir- 
cles. From these plots we can see how the Hen A4686/H/3 
ratio depends on age, metallicity and ionization parameter. 
The lowest metallicity considered in this work is Z = 0.02 
Zq and for this metallicity there is a strong discrepancy 
between the model predictions and the observed data but 
predictions for other metallicities at 4 Myr agree well with 
the observed He II A4686/H/3 ratios. The models can predict 
the Hell A4686 emission line ratio, but only for instanta- 
neous bursts with metallicity of 20% solar and above, and 
only for ages of ~ 4 — 5 Myr, the period when the extreme- 
ultraviolet continuum is dominated by emission from WR 
stars. For burst ages younger than 4 Myr and older than 



6 Myr, and for models with a continuous star formation 
(not shown here), the softer ionizing continuum results in 
an emission spectrum that has too weak Hen lines to be 
consistent with the observational data. 



4.3 The effect of binary evolution on the He II 
4686 emission 

The Stb99 models consider single-star evolution only, but 
it is well-known that massive stars are frequently found in 
binaries and higher order systems which can have a major 
effect on the evolution of massive stars. To explore this pos- 
sibility we compare the observed Hell A4686/H/3 ratio to 
the prediction of the Binary Population and Spectral Syn- 
thesis (BPASS, Eldridge et al 2008, 2009, 2011) model. The 
BPASS code includes a careful treatment of the effect of 
binary evolution on massive short lived stars, and Eldridge 
et al found that including massive binary evolution in the 
stellar population leads to WR stars forming over a wider 
range of ages up to 10 Myr which increases the UV flux at 
later times. In Figure[TT]and Figure[l2]we show the observed 
He IIA4686/H/3 ratio in comparison with their instantaneous 
burst binary and single-star population models, respectively. 
Each panel shows four different metallicities (0.05, 0.2, 0.4,1) 
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Zq . The grey dashed-dotted line shows the N2 PP04 metal- 
licity calibration for 12 + logO/H = 8.2. We should note 
that the lowest metallicity in these plots is a factor of 2.5 
higher than that of Stb99 (0.02 Z Q ). 

We get the highest value for the Hen A4686/H/3 ra- 
tio at 3 Myr and the period with elevated Hen A4686/H/3 
lasts longer. Comparing Figure [TT] and Figure 1121 there is 
not a striking difference in the predicted peak ratio for 
Hell A4686/H/3 between single-star and binary population 
models but clearly the binary model predict an elevated ra- 
tio for a longer period of time than the single-star models. 



5 THE ORIGIN OF NEBULAR HE II 4686 
EMISSION 

The preceding discussion makes it clear that in standard 
models for stellar evolution, only during the phases where 
the ionizing spectrum is dominated by WR stars will we 
see str ong Hen A4686 emission. This was already pointed 
out by ISchaererl (| 19961 ) and our results are in good agree- 
ment with th at work as well as a number of other pre- 
vious works dSchaerer fc Vaccal 1 19981 : [G uscva et all l200d : 
iThuan fc Izotovll2005l l. " 

What has been less studied is the direct test of this 
prediction, namely to ask wheth er WR featur e s are seen 
whenever Hen A4686 is observed. iKehrig et. all (|201lT ) and 
iNeugent fc Massevl (|201lh have recently studied individual 
star forming regions with He n A4686 emission in the Local 
Group and have shown that while most He II -emitting re- 
gions also show evidence of WR stars, not all do. In more dis- 
tant galaxies, previous efforts have primarily looked at He n 
emission in samples selected for other purposes. As an exam- 
ple, B08 looked at Hen galaxies in SDSS DR6 and studied 
whether they showed WR features, but their sample selec- 
tion was not optimised for finding galaxies with Hel l A4686 
emission. The excellent analysis of lThuan fc Izotovl (j2005l ) 
is another example, where the focus was on BCD galaxies. 
Since the present sample is selected purely on the presence 
of He n A4686, we can carry out this study with more relia- 
bility but we need to check for WR features in the spectra 
of our galaxies. 

The presence of WR stars can be recognised via the 
WR bumps around A4650 A(blue bump) and A5808 A(red 
bump). The former is a blend of Hell A4686 and several 
met al lines while the latter is cause d by CivA5801 — 12 
(see ICrowther ARAfcA review! \200H ). B08 and references 
therein). We have therefore inspected the spectra in our star- 
forming sample following the methodology described in B08 
(see Table [3]). Each spectrum is assigned a classification of 
(no WR), 1 (possible WR), 2 (very probable WR), or 3 
(certain WR). We will consider any spectrum with class 1, 
2 or 3 to show WR features. To check the reliability of these 
classifications we can check whether duplicate observations 
of the same object are given the same classification. There 
are four sets of duplicate observations (0417-51821-513 Class 
2, 0418-51817-302 Class 3, 0418-51884-319 Class 2), (0455- 
51909-073 Class 0, 0456-51910-306 Class 0), (0266-51602- 
089 Class 1, 0266-51630-100 Class 0) and (0308-51662-081 
Class 3, 0920-52411-575 Class 3). Only for 0266-51602-089 
and 0266-51630-100 is there uncertainty whether spectrum 
show evidence of WR features or not — we choose to keep 



the original classifications. This is in agreement with, but 
somewhat better than what we find for the full sample of 
WR galaxies from the B08 sample where the RMS classi- 
fication uncertainty from duplicates is 0.4 classes without 
any apparent dependence on the median S/N of the spectra 
down to S/N~ 10; at lower S/N we do not have sufficient 
numbers of duplicate observations to make a statement. In 
total we find that 116 of objects show WR features. 

To check whether the non-detection of WR features is 
due to a low S /N, we have created stacked spectra of galaxies 
with and without WR features. The result of this exercise is 
shown in in Figure [13] The black line is the stack of spectra 
that show WR features while in red we show the stack of 
non-WR galaxies. We see no sign of WR features in the 
stack, further supporting the notion that this class of objects 
show no signs of WR features. 

We now turn to explore whether there are physical dif- 
ferences in the galaxies showing WR features or not. Figure 
1141 shows the ratio of Hell A4686/H/3 versus oxygen abun- 
dance for galaxies with (blue circles) and without WR fea- 
tures (triangles). The symbols for the non-WR galaxies, here 
and in the following, are coloured red for 12 + log O/H < 8.2, 
and orange otherwise. At oxygen abundances lower than 8.2 
we see there is a trend of increasing Hell A4686/H/3 ratio 
towards lower metallicities. We interpret this as being due 
to a harder ionizing continuum at lower metallicities (e.g. 
IThuan fc Izotovll2005l . and references therein). We also see 
a trend of increasing Hell A4686/H/3 towards higher metal- 
licity for 12 + log O/H > 8.7. It is less clear what causes this, 
but since these are more massive systems with higher star 
formation rates and with stronger stellar winds due to their 
higher metallicity, it is likely that what we are seeing is due 
to an increased contribution of shocks and/or a low- level 
AGN contamination. 

Another striking result is that for 12 + log O/H > 8.2, 
essentially all He II -emitting galaxies show WR features, in 
agreement with what one would expect from the models 
discussed in the previous section. At high metallicity there 
appear to be some systems that show a high HellA4686/H/3 
ratios but no sign of WR stars. Since these are more mas- 
sive systems, and fall intermediate between the SF and AGN 
groups in the Hell A4686/H/3 versus [Nil] A6584/Ha dia- 
gram, we interpret this as a likely sign of a low-level AGN 
contribution. 

In Figure [15] we show the fraction of galaxies that show 
WR features as a function of metallicity. To calculate this 
we draw a number of random realisations of the data us- 
ing the uncertainty estimates on the oxygen abundance to 
draw a random realisation. We also draw a random realisa- 
tion of the WR classification assuming a random uncertainty 
in the classification of 0.4 as determined from the analysis 
of duplicate spectra as discussed above. We carry out 101 
random realisations for each of the 101 bootstrap repeti- 
tions and calculate the median fraction of galaxies showing 
WR features in each metallicity bin as well as the 16%-84% 
scatter around the median which is shown by the error bars 
in Figure 1151 What is clear is that there is a transition at 
around an oxygen abundance of 12 + log O/H ~ 8.2 ± 0.1. 
The uncertainty of 0.1 dex encapsulates the fact that the ex- 
act value of this transition abundance depends somewhat on 
the metallicity calibration adopted and 0.1 dex corresponds 
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Figure 13. This plot shows a comparison between two stacked spectra of galaxies with Hen emission. The black spectrum shows the 
result of stacking the spectra of all galaxies in our sample that show WR features. In this case we can clearly see the blue WR bump 
(gray shaded region). The red spectrum shows the result of stacking the spectra of all galaxies that show no WR features. Despite the 
increase in S/N we see no sign of WR features, strengthening our claim that this class of galaxies show no signs of WR features. 



to the scatter found when using the different calibrations 
discussed earlier. 



As we showed in earlier, e.g. Figure 1101 the 
He IIA4686/H/3 ratio depends strongly on the age of the star- 
bursts so it is reasonable to ask whether the systems with- 
out WR features are systematically younger or older than 
the systems that show WR features. In Figure [TB] we test 
this by plotting the Hell A4686/H/3 ratio versus EW(H/3) 
for the sample, where we take the EW(H/3) as a proxy for 
starburst age. We show galaxies with WR features by black 
and blue circles at high and low metallicities, respectively. 
Galaxies without WR features are as before shown by or- 
ange and red triangles at high and low metallicities. We see 
the Hell A4686/H/3 ratios decrease with EW(H/3) for both 
WR and non-WR objects. However, galaxies without WR 
features show higher ratios than galaxies with WR feature 
for the same EW(H/3), especially at lower EW(H/3). Alter- 
natively one might say that at a fixed Hell A4686/H/3 the 
systems without WR features have a higher EW(H/3), or 
with our assumption, a younger starburst. It is not possible 
with our data to disentangle these two possibilities. 
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Figure 14. The Hell A4686/H/3 ratio versus oxygen abundance 
for the sample. The blue solid circles show the location of galaxies 
with WR features. Galaxies which do not show WR features are 
indicated by red and orange triangles at lower and higher oxygen 
abundances than 12 + logO/H = 8.2, respectively. 
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Figure 15. The fraction of objects with detected WR features 
in the He II sample as a function of gas-phase oxygen abundance. 
The points show the median fraction in each abundance bin and 
the error bars the 16%-84% scatter around the median (see text 
for details). While essentially all high-metallicity star- forming 
galaxies with HellA4686 nebular emission show WR features, this 
fraction drops rapidly at metallicities below 12 + logO/H ss 8.2. 
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Figure 16. This plot shows the dependence of the He IIA4686/R3 
ratio on the EW(H/3) and metallicity. We show galaxies with 
WR features by black and blue circles at high and low metal- 
licities, respectively. Galaxies without WR feature are shown by 
orange and red triangles at high and low metallicities. We see 
the Hen A4686/H/3 ratios decrease with EW(H/3) for both WR 
and non-WR objects. However, galaxies without WR features 
show higher ratios than galaxies with WR feature for the same 
EW(HjS), especially at lower EW(H/3)s, alternatively one could 
read this to say that they have a high EW(H/3) (young age) for 
a given Hell A4686/H/3 ratio. 



6 WHY ARE THERE GALAXIES WITH HE II 
4686 EMISSION BUT NO WR FEATURES? 

The models that we discussed previously agree that WR 
stars are the source of the hard ionizing photons necessary 
to produce the Hen A4686 emission. It is therefore a puzzle 
why some of the galaxies with nebular He II A4686 emission 
do not show stellar WR feature s. This lack of WR fea tures 
has been pointed out before (e.g lThuan fc Izotovi 2005), but 



this is the first time a clear trend with metallicity appears. 
In this section we therefore turn to discuss some possible 
reasons for this lack of WR features. 

1. Differences in the S/N of the spectra 

To detect the WR features we need fairly high S/N in the 
continuum as they are broad and weak features. In the top 
panel in Figure [17] we show the relationship between the 
equivalent width of the WR blue bump and the median S/N 
in the continuum for WR galaxies in SDSS DR7. We see that 
a S/N of 10 is sufficient to detect even very weak features. 
The distribution of the S/N for WR and non-WR objects 
in the sample versus their metallicities in the bottom panel 
shows that 20 out of 70 of the non-WR galaxies at 12 + 
log O/H < 8.2 have S/N less than 10. Since the total number 
of galaxies at 12 + log O/H < 8.2 is 115, if all the 20 low 
S/N galaxies are assumed to have WR features this brings 
the fraction of galaxies with WR features in this metallicity 
range from 30% to 43%. Thus low S/N could cause us to 
underestimate the WR fraction by at most ~ 15%. But the 
fact that the co-added spectrum in Figure [T3] which has a 
higher S/N shows no WR features is suggestive that S/N 
might not be the problem for non-detection of WR features. 

Furthermore even adding 15% would still mean that less 
than half He II emitters at low metallicity would show WR 
features. We also studied the difference in WR classification 
of duplicate observations as a function S/N and as remarked 
earlier, we saw no trend with S/N for more uncertain clas- 
sification down to a S/N of 10. The dataset is insufficient 
to test this at lower S/N. Thus our conclusion is that it is 
unlikely that the systematic absence of WR features in low 
metallicity objects is due to low S/N in the spectra. 

2. Weak lined WN stars 

One possible reason for the lack of WR features in the 
most metal poor galaxies, is that the WR lines are too 
weak to be seen. It is we ll known (e.g. IConti et al.l 1 19891 ; 
ICrowther fc Hadfieldl2006l ) that WR stars in the SMC have 
narrower and le ss luminous lines tha n equi valent stars in 
the Milky Way. ICrowther fc Hadfieldl l|2006h . for instance, 
find that He II lines in WR stars in the SMC are typically 
a factor of 4-5 weaker than the Milky Way and hence one 
would expect that in galaxies at the same distance, it would 
be harder to see WR features in lower metallicity systems. 
This is however countered by the fact that low metallicity 
systems on average are closer, and hence the SDSS fibre 
subtends a smaller physical size. Since WR emitting regions 
typically are small, they do not fill the 3" aperture in more 
distant galaxies, which means that the contrast of the WR 
features is being enhanced in low redshift systems. This is 
reflected in the fact that in the WR survey by B08, the 
equivalent width of the WR features in low metallicity sys- 
tems (12 + Log O/H < 8.25) is slightly higher than that 
in metal rich objects (12 + Log O/H > 8.5), a mean of 5.1 
A vs 4.1 A. Thus the increased contrast appears to approx- 
imately cancel out the decrease in line luminosity leading 
to a fairly constant detection potential with redshift. Thus 
we do not believe that this is the cause of the dearth of 
WR features in the low metallicity Hell emitting galaxies, 
at least down to 12 + log O/H ~ 8. At very low metallici- 
ties, 12 + log O/H ~ 7.5, Eldridge & Stanway (2009) found 
that in their models the WR features become very weak 
and if those results are correct it should be extremely hard 
to detect WR features in those galaxies. We do however 
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caution that I Zw 1 8 has very prominent WR features (e.g. 
iLegrand etafll 19971 ). so at least some extremely metal poor 
galaxies do show clear WR features. Furthermore, even if 
we are unable to detect WR features at the very lowest 
metallicities, the same models predict strong WR features 
at 12 + logO/H > 8, thus this is not a sufficient explanation 
for the result in Figure [151 

3. Shocks 

iThuan fc Izotovl £2005) studied the hard ionizing radiation 
in very metal poor BCD galaxies in the local Universe and 
concluded that fast radiative shocks could be responsible for 
the nebular Hen A4686 emission. Therefore, another pos- 
sibility is that there is a contribution to the Hel l A4686 
from shocks in the ISM (|Dopita fc Sutherland! Fl996l ). That 
He n A4686 can have some contribution from shocks is plau- 
sible, but whether it can explain the systematic lack of WR 
features at low metallicity is less clear. 

One question is whether shock models can reproduce 
the line luminosity i n Hell at low me tallic ities. Using 
the pr edictions of the iDopita et all (|2005l ) and [Ah en et al.l 
(2008) shock models for the Hell A4686/H,5 ratio versus 
[Nil] A6584/Ha, we found that we can only obtain a ratio 
comparable to the observed one for objects at high metal- 
licity. A second issue centers on the observation that shocks 
would most likely come from supernovae and stellar winds, 
but the latter are though to be weaker at low metallicities. 
Shocks can also be induced by outflows from starbursts (e.g. 
galactic winds) and mergers but only two galaxies in our 
sample are interacting and show a perturbed morphology 
and we see no clear difference between the low- and high- 
metallicity subsamples. 

4. X-ray binaries 

Another candidate source for He II ionization that has been 
discussed in previous studies of He H A4686 emittin g nebu- 
lae are massive X-ray binaries |Garnett et al. 199 ll ). While 
these are likely present in many active star forming regions, 
the question here is why massive X-ray binaries should be 
more common at low metallicity. If the He + ionizing pho- 
tons at low metallicities come from X-ray binaries, we would 
expect an increasing binary fraction with decreasing metal- 
licity. Without a theoretical justification for this, we consider 
an increased abundance of X-ray binaries at low metallicity 
to be an unlikely explanation for the trend seen in Figure [151 

5. post-AGB star s 

iBinette et alj l| 19941 ) demonstrated that photoionization by 
post-AGB stars can produce nebular Hell A4686 emission. 
After about a few 10 7 years (i.e, after massive stars disap- 
pear), the ionizing radiation comes from post-AGB stars. 
So, the ionization from post-AGB stars become more im- 
portant in more evolved systems and this is not the case 
for our objects especially not for objects at low metallicities 
(c.f. Figure [To)) . 
6. Spatial offset 

One possible explanation for non detection of WR features is 
that there could be a significant spatial separ ation between 
the W R stars and the region emitting He II . iKehrig et. all 
(2008) saw indeed such a spatial separation based on inte- 
gral field spectra of II Zw 70. The found that the location 
of the WR stars and the He II A4686 emission appear to be 
separated by ~ 80pc. 

Similarly, llzotov et al.l (|2006bl ) studied two-dimensional 
spectra of an extremely metal-deficient BCD galaxy SBS 
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Figure 17. Top panel: The relationship between the equivalent 
width of the WR blue bump and the median S/N in the contin- 
uum for WR galaxies in SDSS DR7. We see that a S/N of 10 is 
sufficient to detect even very weak features. Bottom panel: The 
distribution of the median S /N in the continuum for WR and non- 
WR objects in the sample versus their oxygen abundance, we see 
that 20 out of 70 of the non-WR galaxies at 12 + logO/H < 8.2 
have S/N less than 10. See text for more discussion. 



0335-052E and showed that the HellA4686 emission line was 
also offset from the near evolved star clusters but in their 
case, by studying the kinematical properties of the ionized 
gas from the different emission lines they suggested that the 
hard ionizing radiation responsible for the He II A4686 emis- 
sion was not related to the most massive youngest stars, but 
rather was related to fast radiative shocks. 

If this offset between the WR stars and the region emit- 
ting He II might be an explanation for our non-detections of 
WR features, it would mean that such a spatial separation 
is much more common at low metallicity which is rather sur- 
prising, since stellar winds are thought to be considerably 
weaker there, however low metallicity galaxies are also on 
average closer so the SDSS fibre subtends a smaller physical 
scale so a smaller volume would need to be blown out by 
the wind (see Figure 18). To get a more quantitative esti- 
mate we calculate the gravitational binding energy of a cloud 
with radius 1.5" (SDSS aperture radius) at the redshift of 
each non-WR object. We assume a hydrogen density of ~ 
50 cm -3 . The median energy require to excavate a hole of 
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this size, is of the order of 10 55 erg for our full sample. For 
the lowest redshift objects the energy requirement is a much 
more manageable 10 49 erg and thus in these cases the ab- 
sence of WR features in the spectrum could be due to a spa- 
tial offset from the He II emitting region. At higher redshift 
the energetics makes this a much less likely explanation. The 
possibility does however warrant further examination and to 
test this we are undertaking a spectroscopic follow-up of a 
subsample of these sources. 

7. Chemically homogeneous stellar evolution 

A final explanation could be that the stellar populations at 
very low metallicities can have much higher temperatures 
than is currently expected in models. This would be the case 
if some stars ro t ated fast enough to evolv e homogeneously 
ilMaederl 1 19871; iMevnet fc Maederl 120071: lYoon fe Langerl 
120051 ; lYoon et al.ll2006l ; ICantiello et al.ll2007h . In that case 



we can get a higher continuum at 228 A and correspond- 
ingly a higher HellA4686/H/3 ratio in comparison with non- 
homogeneous stellar evolution models. An appealing aspect 
of this, speculative, explanation is that homogeneous evo- 
lution is predicted to be more common at low metallic- 
ity. There are however currently no studies of the nebu- 
lar He n A4686 line in the literature. lEldridge fc Stanwavl 
(|2011bT ) looked at the effect of quasi-homogeneous evolu- 
tion in their binary models and showed that including it led 
to strengthened WR features at low redshift, in contrast to 
what we need. However that still leaves the possibility open 
that there is a period where strong nebular He 11 A4686 is 
seen but no WR features as that has not yet been tested. 



7 CONCLUSION 

We have presented a sample of rare star-forming galaxies 
with strong nebular Hen A4686 emission spanning a wide 
range in metallicity. We have derived physical parameters 
for these galaxies and showed that emission line models that 
can reproduce the strong lines in the galaxy spectra are not 
able to predict the observed ratio of Hen A4686/H/3 at low 
metallicities. In agreement with previous studies we found 
that current models for single massive stars are able to re- 
produce the He 11 A4686/H/3 ratio in galaxies in our sample, 
but only for instantaneous bursts of 20% solar metallicity or 
higher, and only for ages of ~ 4 — 5 Myr, the period when 
the extreme- ultraviolet continuum is dominated by emission 
from WR stars. For stars younger than 4 Myr or older than 
5 Myr, and for models with a constant star-formation rate, 
the softer ionizing continuum results in Hell A4686/H,9 ra- 
tios typically too low to explain our data. Including massive 
binary evolution in the stellar population analysis leads to 
WR stars occurring over a wider range in age which leads 
to acceptable agreement with the data at all metallicities 
sampled as long as WR stars are present. 

However, the most notable result of our studies is that 
a large fraction of the galaxies in our sample do not show 
WR features and this fraction increases systematically with 
decreasing metallicity. We find that 70% of galaxies at oxy- 
gen abundances lower than 8.2 do not show WR features in 
their spectra. We discussed a range of different mechanism 
responsible for producing He 11 A4686 line apart from WR 
stars in these galaxies and conclude that spatial separation 
between WR stars and the region emitting Hell emission 



can be a possible explanation for non-detection of WR fea- 
tures in these galaxies. Moreover, if the stellar population 
models at very low metallicities can have much higher tem- 
peratures than is currently expected in models, as would 
for instance be the case if some stars rotate fast enough to 
evolve homogeneously, then such models might explain the 
origin of the He 11 A4686 line and also the metallicity trend of 
the Hell sample better. We will explore these possibilities 
in a future paper. 
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APPENDIX A: FITTING MODELS TO THE 
EMISSION LINES 

The CL01 model grid is calculated by varying the model 
parameters, U, ionization parameter at the edge of the 
Stromgren sphere, tv, total l/-band optical depth, £, the 
dust-to-metal ratio of ionized gas, /x, the fraction of the total 
optical depth in the neutral ISM contributed by the ambi- 
ent ISM, and Z, the metallicity, over a certain range for 221 
unequally spaced time steps from t = to t = 20Gyr (see 
Table |J). 

In this paper we use the interpolated model grid of various 
luminosities for 50 time-steps from B04. This includes a to- 
tal of 2 x 10° different models. To fit to the data we adopt 
the Bayesian methodology described by Ka03. We obtain 
the PDF of every parameter of interest by marginalisation 
over all other parameters. The resulting PDF is used to es- 
timate confidence intervals for each estimated physical pa- 
rameter. We need to fit to at least five strong emission lines, 
[On] A3727,3729, H/3, [Om] A5007, Ha, [Nil] A6584 to get 
a good constraint on the parameters. We take the median 
value of each parameter to be the best estimate of a given 
parameter. 

In Figure IAT1 we illustrate our technique by showing the 
effect of adding lines on the PDFs of parameters when we 
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Figure 18. Images of nine galaxies which do not show WR features in their spectra while have strong Hen A4686 emission. The sizes 
of boxes are 50 " X 50 ". In these images north is up, east to the left. Red circles show the size of SDSS fibre aperture (3") and yellow 
boxes show 100 pc around the center of the SDSS fibre. 

fit a model to the data. We start with [O n] A3727, 3729 and 
show how we get more well defined PDFs for the indicated 
parameters as we add the emission lines indicated on the left. 
We show the PDFs for dust attenuation parameter in V - 
band, gas phase oxygen abundance, ionization parameter, 
dust-to-metal ratio of ionized gas and the conversion factor 
from Ha and [O n] luminosity to star formation rate (see 
CL01 for further details), for one object in our sample. The 
dust-to-metal ratio, £, is hard to constrain, except at high 
metallicity. 
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t v 12 + Log 0/H Log U f rj[Ha] ry[0 II] 




1 2 3 7.5 8.0 8.5 9.0 -4.0 -3.5 -3.0 -2.5 0.1 0.2 0.3 0.4 5.5 6.0 6.5 7.0 7.55.5 6.0 6.5 7.0 7.5 



Figure Al. We show how the PDFs for dust attenuation parameter in V-band, gas phase oxygen abundance, ionization parameter 
of ionized gas, dust-to-metal ratio and Ha and [Oil] efficiency factors change when we add more emission lines to the fit. £ is hard to 
constrain, except at high mctallicity. 
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Table 2: The positions and the identifications of galaxies with neb- 
ular Hen A4686 emission. OKSN is for galaxies having SN > 3 in 
their [Om] A5007, H/3, [Nil] A6584 and Ha lines. See B04 for the 
BPT classification. The full table is available in electronic form in 
http: / /www.strw. leidenuniv.nl/~shirazi/SB 01 1 / 1 
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Table 5: The positions and identifications of the sample of star-forming 
galaxies with nebular He II A4686 emission. 



Plate-MJD-FiberlD 


Z 


a(J2000) 


5 


ID 


log Hen /H/3 


log [Nil] /Ha 


12 + logO/H 


Features 


Class 


Other Names 


0752-52251-340 


0.00291 


+00 


09 


53.09 


+15 44 04.80 


51 


-1.35 


-1.67 


7.92+0.09 


Non-WR 





N/A 


0390-51900-291 


0.09437 


+00 


17 


28.29 


-00 56 24.98 


13 


-1.57 


-0.56 


9.27+0.05 


WR 


2 


N/A 


0390-51900-445 


0.09840 


+00 


21 


01.03 


+00 52 48.08 


14 


-1.96 


-1.10 


8.83+0.04 


WR 


2 


UM 228 


0753-52233-094 


0.01424 


+00 


24 


25.95 


+14 04 10.65 


52 


-1.98 


-1.36 


8.10+0.12 


WR 


3 


N/A 


0418-51884-319 


0.01791 


+00 


32 


18.59 


+15 00 14.17 


19 


-1.97 


-1.52 


8.01+0.07 


WR 


2 


SHOC 022 


0691-52199-389 


0.19554 


+00 


42 


52.31 


+00 27 30.09 


48 


-0.75 


-0.82 


8.88+0.05 


WR 


2 


N/A 


1905-53706-628 


0.15928 


+00 


49 


13.88 


+00 24 01.99 


138 


-1.64 


-2.05 


7.85+0.06 


Non-WR 





N/A 


0394-51913-075 


0.16742 


+00 


55 


27.46 


-00 21 48.77 


15 


-1.76 


-0.89 


8.90+0.04 


WR 


2 


N/A 


0695-52202-137 


0.00556 


+01 


15 


33.82 


-00 51 31.17 


49 


-2.32 


-1.27 


8.31+0.35 


WR 


3 


NGC 0450 


2329-53725-536 


0.00975 


+01 


25 


34.19 


+07 59 24.40 


169 


-2.02 


-2.22 


7.80+0.04 


Non-WR 





N/A 


0429-51820-495 


0.05662 


+01 


47 


07.04 


+13 56 29.29 


20 


-1.95 


-1.52 


8.02+0.05 


Non-WR 





N/A 


0666-52149-331 


0.15155 


+01 


59 


53.07 


-08 13 48.99 


44 


-1.73 


-1.43 


8.04+0.06 


Non-WR 





SHOC 099 


1073-52649-409 


0.03994 


+02 


13 


06.62 


+00 56 12.44 


73 


-1.97 


-1.58 


8.09+0.05 


Non-WR 





UM 411 


0667-52163-634 


0.00492 


+02 


15 


13.98 


-08 46 24.39 


45 


-1.12 


-1.81 


7.74+0.09 


Non-WR 





SHOC 111 


0456-51910-306 


0.08222 


+02 


40 


52.20 


-08 28 27.43 


24 


-1.60 


-1.63 


8.02+0.04 


Non-WR 





SHOC 133 


1070-52591-072 


0.00421 


+02 


42 


39.86 


-00 00 58.64 


72 


-1.41 


-0.33 


8.93+0.25 


Non-WR 





M077 


0456-51910-076 


0.00456 


+02 


48 


15.94 


-08 17 16.51 


23 


-2.11 


-1.90 


7.87+0.08 


Non-WR 





N/A 


1666-52991-310 


0.03992 


+03 


14 


31.71 


+41 05 25.95 


117 


-1.95 


-1.40 


8.08+0.19 


WR 


1 


N/A 


0413-51821-480 


0.06868 


+03 


17 


43.12 


+00 19 36.84 


16 


-0.77 


-0.54 


9.07+0.11 


WR 


2 


N/A 


1733-53047-528 


0.07881 


+07 


29 


30.29 


+39 49 41.62 


121 


-1.97 


-0.93 


8.90+0.04 


Non-WR 





N/A 


1922-53315-588 


0.06978 


+08 


06 


19.50 


+19 49 27.31 


139 


-2.05 


-1.19 


8.83+0.04 


WR 


2 


N/A 


0761-52266-361 


0.04592 


+08 


22 


27.44 


+42 23 31.14 


53 


-1.09 


-0.37 


9.25+0.07 


WR 


1 


N/A 


1267-52932-384 


0.04722 


+08 


23 


54.97 


+28 06 21.75 


78 


-2.16 


-0.95 


8.88+0.04 


WR 


1 


N/A 


0548-51986-503 


0.00706 


+08 


26 


04.80 


+45 58 07.36 


30 


-2.07 


-0.76 


8.64+0.26 


WR 


1 


UGC 04393 


2425-54139-372 


0.01974 


+08 


29 


32.66 


+14 27 06.92 


173 


-2.00 


-1.35 


8.13+0.20 


Non-WR 





N/A 


0445-51873-404 


0.00245 


+08 


37 


43.48 


+51 38 30.26 


21 


-2.38 


-1.71 


7.95+0.08 


Non-WR 





MRK 0094 


0828-52317-148 


0.14746 


+08 


38 


43.64 


+38 53 50.50 


56 


-1.88 


-1.32 


8.69+0.05 


WR 


1 


N/A 


2278-53711-411 


0.07219 


+08 


40 


00.37 


+18 05 31.01 


168 


-2.02 


-1.53 


7.98+0.07 


Non-WR 





N/A 



g3 

ix> 

. 

s 



S5 
S 
S 



Table 5: continued. 
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